clear all; close all; clc;
Horb=1.91E4;
T=60; %шаг модели
Rz=6.4E3;
Rorb=Horb+Rz;
fprintf ('Rorb=%f km\n',Rorb);
torb=0:T:6*3600;
tob=12*3600;
alf=2*pi*torb/tob;
Rsv_x=Rorb*cos(alf);
Rsv_y=Rorb*sin(alf);
fprintf ('Rsv_x=%f км\n', Rsv_x);
Rr=[0,Rz];
Rv_x=Rsv_x-Rr(1);
Rv_y=Rsv_y-Rr(2);
V_x=diff(Rv_x)/T;
V_y=diff(Rv_y)/T;
figure(1);
V_mod=sqrt(V_x.^2+V_y.^2); %Модуль вектора скорости
plot(torb(2:end),V_mod);
R_mod=sqrt(Rv_x.^2+Rv_y.^2); %Модуль радиус вектора
Rn_x=Rv_x./R_mod;
Rn_y=Rv_y./R_mod;
Vv=Rn_x(2:end).*V_x+Rn_y(2:end).*V_y; %Скалярное произведение вектора скорости на единичный вектор по линии визирования - проекция вектора скорости на линию визирования
f0=1.6e9;
c=2.99e5;
fd=f0*Vv/c;
k1=0; k2=0;
for k=2:length(torb)
    if Rv_y(k)>0
        k1=k1+1;
        Dop(k1)=fd(k-1);
        x1(k1)=Rv_x(k);
        y1(k1)=Rv_y(k);
    else
        k2=k2+1;
        x2(k2)=Rv_x(k);
        y2(k2)=Rv_y(k);
    end
end
figure(2);
hist(Dop);
